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ABSTRACT 

We conducted a deep survey of resolved sub-systems among wide binaries with solar-type 
components within 67 pc from the Sun. Images of 61 stars in the K and H bands were obtained 
with the NICI adaptive-optics instrument on the 8-m Gemini-South telescope. Our maximum 
detectable magnitude difference is about 5™ and 7.8™ at 0.15" and 0.9" separations, respectively. 
This enables a complete census of sub-systems with stellar companions in the projected separation 
range from 5 to 100 AU. Out of 7 such companions found in our sample, only one was known 
previously. We determine that the fraction of sub-systems with projected separations above 5 
AU is 0.12 ± 0.04 and that the distribution of their mass ratio is flat, with a power-law index 
0.2 ± 0.5. Comparing this with the properties of closer spectroscopic sub-systems (separations 
below 1 AU), it appears that the mass-ratio distribution does not depend on the separation. 
The frequency of sub-systems in the separation ranges below 1 AU and between 5 and 100 AU is 
similar, about 0.15. Unbiased statistics of multiplicity higher than two, advanced by this work, 
provide constraints on star-formation theory. 



Subject headings: stars: binaries 

1. Introduction 

The quest for exo-planets has attracted atten- 
tion to stars in the solar neighbourhood. A com- 
plete understanding of their formation requires a 
quantitative and predictive description of stellar 
multiplicity. Empirical data on stellar multiplicity 
serve as a benchmark for star formation theories. 



for example to compare wit h recent larg e-scale 
hydrodynamical simulations ( Bate I 120091 ). Al- 



though multiplicity statistics of solar-type dwarfs 
were assumed to be known since the survey of 



Duquennov fc Mavor I (|l99lh . in fact that survey 



did not constrain multiplicities higher than two 
because (i) the sample of 164 stars was too small 
and (ii) the discovery of wide binaries re l ied on 



Based on observations obtained at the Gemini Obser- 
vatory (Program ID GS-2009B-Q-49), which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership: the National Sci- 
ence Foundation (United States), the Science and Tech- 
nology Facilities Council (United Kingdom), the National 
Research Council (Canada), CONICYT (Chile), the Aus- 
tralian Research Council (Australia) , Ministerio da Ciencia 
e Tecnologia (Brazil) and Ministerio de Ciencia, Tecnologia 
e Innovacion Productiva (Argentina). 
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shows that iDuquennov fc Mavor I ( |l99l[ ) missed 



traditional visual techniques. lRa.ghavan 



half of higher-order multiples within 25 pc from 
the Sun. 

Several studies of large samples of nearby stars 
are now in progress. Radial-velocity (RV) tech- 
niques ensure a secure detection of all short-period 

stella r companions (see summary in lGrether fc Lineweaver 



20061 ). However, for wider separations (> 10 AU) 
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RV variations are too slow or too small, and such 
companions are better discovered by direct imag- 
ing. Sev eral recent deep imaging su r veys, summa- 
rized by iMetchev fc Hillenbrandtl (|2009l) . aimed 
to detect substellar or planetary-mass companions 
to nearby stars, obtaining statistics of stellar com- 
panions as a by-product. The typical sample size 
of imaging surveys is still small, on the order of 
100 targets. 

Most imaging and RV surveys for low-mass 
compani ons explicitly a.void known visual bina- 
ries (e.g. Chauvin et al. llioiol ). so important star- 
formation effects may not be represented in their 
samples. For example, catalogued triples show 
that the mass-ratio distributions in close inner 
sub-systems discovered spectroscopically and in 
wider "visual" sub-systems are radically different 
(see Section 5). One of our goals is to investigate 
if this is a selection effect or a genuine feature of 
the star-formation process. 

In order to improve the multiplicity statistics 
of a well-defined and sufficiently large distance- 
limited sample, we focus on ^5000 dwarfs with 
0.5 < y - / < 0.8 within 67 pc from the Sun. 
According to preliminary estimates, there should 
be some 400 triples and 100 quadruples in this 
sample. Surveying such a large sample to de- 
tect those multiples with a reasonably deep and 
known completeness limit is a large task; here we 
address only a small part of it. Instead of avoid- 
ing known binaries, we purposefuly selected wide 
visual binaries from this sample. Some of these 
targets are also monitored in RV, giving comple- 
menta ry constraints on sub- systems with short pe- 



riods (jPesidera et al. 1 120061 ). 



We imaged the chosen wide binaries with Adap- 
tive Optics (AO) to study the properties of in- 
ner sub-systems. AO imaging on 8-m class tele- 
scopes provides nearly diffraction-limited resolu- 
tion of 0.06" in the near-infrared, which corre- 
sponds to a projected orbital radius of 4 AU for the 
farthest target stars in our sample, filling an im- 
portant gap between RV and seeing-limited imag- 
ing surveys. 

Our sample is detailed in Section [51 Observa- 
tions, data reduction, and detection limits are de- 
scribed in Section |3l the results are presented in 
Section 131 Sections [5l and [6l contain the discussion 
and conclusions, respectively. 



2. The sample of wide binaries 

The wide binaries chosen for our survey be- 
long to a large sample of nearby solar-type dwarfs 
[Nsample) sele cted from the latest version of Hip- 
parcos catalog ( van Leeuwen Il2007t) by the follow- 
ing criteria: 

1. trigonometric parallax tthip < 15 mas (within 
67 pc from the Sun, distance modulus < 
4.12™); 

2. color 0.5 < V — I < 0.8 (this corresponds 
approximately to spectral types from F5V 
to KOV); 

3. unevolved, satisfying the condition Afnp > 
9{V — I) — 3.5, where Mnp is the absolute 
magnitude in the Hipparcos band calculated 

with TTHIp. 

There are 5040 stars (i.e. Hipparcos catalog en- 
tries) satisfying these conditions. Some of these 
stars with erroneous parallaxes or colors will be 
eventually removed from the sample, but this con- 
tamination is expected to be minor. The Nsam- 
ple defined above is reasonably complete because 
the cumulative object count versus distance d fol- 
lows the cfi law, deviating from it by only 10% at 
d — 67 pc. 

Visual binaries with separations 5" < p < 
20" and one or both components belonging to 
the Nsample were selected for AO observations 
from the Washington Double-Star Catalog, WDS 



([Mason et al. I [20011 ). Additionally, we require 



each pair to have at least 3 observations and mag- 
nitudes of both components listed in the WDS, 
and that apparent motion of the wide pair be com- 
patible with TTHip. As most targets have large 
proper motions, relative stability of the binary po- 
sition over time permits rejection of optical pairs. 
An additional check was made by placing the com- 
ponents on the {My, V — K) color-magnitude di- 
agram (CMD). There are ~200 wide binaries in 
this list, 101 with negative declinations. 

In Table[Tlwe list basic data on the components 
of the 33 wide pairs observed in this run. Column 
1 gives the Hipparcos numbers of the primary and, 
in some cases, secondary components. The next 
three columns contain tthip, binary separation p, 
and projected separation r = p/tthip. Then fol- 
low the magnitudes of the primary (A) and sec- 
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Table 1 
Data on observed targets 



HIP (A/B) 


WHIP 

mas 


P 


r 

AU 


Component A 
V H 


Component B 

V H K, 


Date 
2000+ 


Rem 


2028/29 


24. 


,0 


16.5 


687 


8, 


.39 


7, 


,46 


7, 


,41 


9.56 


8, 


,30 


8, 


,27 


09, 


,9400 


B 


10579 


17, 


,2 


6.7 


390 


9, 


.43 


8, 


,06 


7, 


,97 


9.97 


8, 


,34 


8, 


,24 


09, 


,9400 




20552 


36, 


,1 


5.6 


155 


6 


.87 


5, 


,37 


5, 


,32 


7.23 


5, 


,64 


5, 


.55 


10, 


,0851 




20610/12 


16, 


,7 


10.5 


629 


8, 


.06 


6, 


,89 


6, 


,84 


8.27 


7, 


,01 


6 


.98 


10, 


,0851 


B 


21963 


18, 


,2 


8.2 


450 


8, 


.13 


6, 


,69 


6, 


,62 


13.30 


8, 


,99 


8, 


,75 


10, 


,0852 




22531/34 


26, 


,0 


12.9 


496 


5 


,61 


4, 


,88 


4, 


,80 


6.24 


5, 


,28 


5, 


,19 


09, 


,9404 


A,CRN 


23926/23 


18, 


7 


10.1 


540 


6 


,83 


5, 


,20 


5, 


10 


10.29 


8, 


,00 


7, 


,90 


10, 


,0851 




24711/12 


15, 


,3 


13.3 


867 


8. 


,46 


7, 


,00 


6, 


,95 


10.57 


8, 


,29 


8, 


,24 


10, 


,0852 




27922 


42, 


,4 


10.4 


246 


7. 


,60 


5, 


,84 


5, 


76 


10.57 


7, 


,40 


7, 


,22 


10, 


,0852 




28790 


37, 


,2 


5.6 


150 


6, 


.02 


4, 


,92 


4, 


,75 


8.98 


6, 


,09 


6, 


,03 


10, 


,0853 




30158 


17, 


,8 


7.0 


390 


8, 


.57 


7, 


,02 


6, 


,91 


10.58 


8, 


,09 


8, 


,03 


10, 


,0853 




32644 


21, 


,2 


4.9 


230 


7, 


.37 


6, 


,05 


5, 


,97 


8.66 


6, 


,61 


6 


,52 


10, 


,0852 




36165/60 


31, 


,2 


17.7 


566 


7, 


.08 


5, 


,92 


5, 


,80 


8.06 


6, 


,75 


6 


.60 


10, 


,0854 




37332/33 


15, 


,1 


15.7 


1043 


7, 


,89 


6, 


,47 


6, 


,41 


10.01 


8, 


,35 


8, 


,25 


10, 


,0855 




37335 


21, 


,8 


6.1 


282 


8 


,93 


7, 


,05 


6, 


,99 


12.38 


8, 


,59 


8, 


,40 


10, 


,0855 




39409 


15, 


,7 


5.1 


328 


9 


,26 


7, 


,57 


7, 


,44 


9.36 


7, 


,56 


7, 


,43 


10, 


,0855 




43947 


26, 


,6 


10.2 


383 


8. 


,44 


7, 


,28 


7, 


19 


9.95 


8, 


,27 


8. 


,20 


10, 


,0857 




44584/85 


17, 


,0 


10,5 


617 


8, 


,09 


6, 


,99 


6, 


,90 


9.19 


7, 


,72 


7, 


,70 


10, 


,0857 


B 


44804 


15, 


,5 


7.4 


478 


8, 


,80 


7, 


,20 


7, 


15 


10.70 


8, 


,45 


8, 


,29 


10, 


,0857 




45734 


17, 


,6 


9.0 


513 


8, 


.41 


6, 


,85 


6, 


,78 


9.66 


7, 


,56 


7, 


,44 


10, 


,0856 




45940 


28, 


,6 


6.8 


237 


8, 


.17 


6, 


,55 


6, 


,43 


11.86 


8, 


,35 


8, 


,07 


10, 


,0856 




46236 


21, 


,4 


19.3 


900 


7, 


.03 


5, 


,77 


5, 


,69 


10.18 


7, 


,82 


7, 


,68 


10, 


,0856 




47839/36 


40, 


,2 


18.8 


468 


8, 


.08 


6, 


,79 


6, 


,70 


8.23 


6, 


,99 


6 


.85 


10, 


,0858 




49520 


16, 


,9 


9.6 


565 


8. 


82 


7, 


,35 


7, 


,24 


8.99 


7, 


,40 


7, 


,30 


10, 


,0858 




50638/36 


83, 


,3 


17.1 


205 


7, 


,67 


5, 


,18 


4, 


,98 


9.67 


8, 


,46 


8. 


,40 


10, 


,0858 




50883 


16, 


,2 


6.7 


413 


7, 


,90 


6, 


,71 


6, 


,62 


13.00 


9, 


,05 


8, 


,98 


10, 


,0858 




55288 


20, 


,8 


9.5 


456 


7, 


,10 


5, 


,91 


5, 


,82 


7.91 


6, 


,32 


6, 


,22 


10, 


,0859 




58240/41 


21, 


,1 


18.9 


895 


7, 


,67 


6, 


,24 


6, 


,13 


7.83 


6, 


,29 


6, 


.24 


10, 


,0859 




59021 


19, 


,4 


6.0 


308 


6, 


.69 


5, 


,24 


5, 


,10 


8.84 


6, 


,71 


6, 


,55 


10, 


,0859 


CRN 


61595 


17, 


,8 


8.4 


472 


8, 


.28 


6, 


,80 


6, 


,68 


10.48 


8, 


,29 


8, 


,12 


10, 


,0860 


A 


64498 


18, 


,0 


9.3 


516 


7. 


,74 


6, 


,40 


6, 


,25 


10.10 


7, 


,98 


7 


,85 


10, 


,0861 




65176 


17, 


,8 


5.0 


282 


7. 


,85 


6, 


,47 


6, 


,40 


8.59 


6, 


,82 


6, 


,68 


10, 


,0861 




67408 


33, 


,9 


11.6 


342 


6, 


.62 


5, 


,29 


5, 


,22 


10.21 


7, 


,14 


7, 


,02 


10, 


,0860 


A, CRN 
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ondary (B) components in the V, H, and K,^ bands 



as lis ted in the WDS and 2MASS (|Cutri et al 
120031 ). The last two columns contain the date 
of observations and remarks where we indicate 5 
cases when only one component was observed and 
3 cases when the coronagraphic mask was used. 

3. Observations and data reduction 
3.1. Observing procedure 

The Near-Infrared Coronagraphic Imager, NICI, 
on the Gemini South telescope is an 85-element 
curva ture AO instrument based on natural guide 



stars (|Toomev fc Ftaclas l2003HChun et a"ni2008l ) 
We used NICI in normal (non-coronagraphic) 
mode, as a classical AO system. However, si- 
multaneous imaging in the two infrared bands H 
and K offered by NICI helps to distinguish faint 
companions from static speckles (the radial dis- 
tance of a static speckle scales with wavelength 
while the position of a real detection does not 
change). The two detectors have 1024^ pixels of 
18 mas (milliarcseconds) size, covering a square 
field of 18". 

To avoid saturation, we used narrow-band 
filters with central wavelengths 2.272 /im and 
1.587 /im in the red and blue channels, respec- 
tively. The minimum possible exposure time is 
0.38 s, still causing saturation for bright targets. 
We observed three of the brightest stars with the 
coronagraphic mask of 0.32" radius, using in this 
instance broadband K and H filters. For the 
remaining bright targets we allowed moderate 
saturation. As the tradeoff between slower ob- 
servations with more complex data reduction in 
coronagraphic mode and slightly saturated PSFs 
is not obvious, we tried both. 

Typically, each target was observed at 5-6 
dithered positions for a total accumulation time 
of about 5 min. The dither offsets were designed 
to take images of both wide-binary components 
in one frame whenever possible {p < 10") while 
closing the AO loop on the primary. Components 
of wider pairs were observed separately, except 
secondaries with V > 12™, too faint to serve as 
AO guide stars. 

Of the two nights allocated to this program, 
one was lost due to a technical problem. Three 
targets were observed on 2009 December 10 in ser- 
vice mode. The remaining objects were observed 



on the night of 2010 January 31 to February 1 in 
classical mode. The seeing during that night was 
good. As measured from the NICI AO loop data, 
it ranged between 0.34" and 0.94", with a median 
of 0.50" (at 500 nm, not corrected to the zenith). 
High AO performance was achieved (see next Sec- 
tion). The transparency was variable (light cir- 
rus). Out of 202 binary components accessible to 
NICI, we were able to observe 61, or 30%. The effi- 
ciency was high, with an average of 12 min per tar- 
get including telescope slew, tuning of the primary 
mirror, acquisition of the object on the science 
camera, closing the AO loop, and data-taking. 

3.2. Data processing 

The data were processed in a standard way us- 
ing custom IDL programs. Individual images writ- 
ten to separate FITS files were combined into data 
cubes, separately in red and blue channels. The 
sky image was obtained for each target as a me- 
dian of the cube, then subtracted from each plane 
of the cube. After correcting for bad pixels and di- 
viding by the fiat field, the images were re-centered 
by cross-correlation and median-combined. 

Most frames contain just one or two bright 
point sources corresponding to the components of 
the wide binary. The position of each point source 
was recorded by "clicking" on its image, then re- 
fined by fitting a 2-dimensional Gaussian to the 
central part of the Point Spread Function (PSF) 
within 5-pixel radius. The ratio of the central in- 
tensity of this Gaussian to the total flux in a 0.9" 
diameter circle gives an estimate of the Strehl ra- 
tio. The typical StrehlfH are 43% and 22% in the 
red and blue channel, respectively, with maximum 
values reaching 54% and 38%. The median Full 
Width at Half Maximum (FWHM) of the PSF 
core was 65 mas and 57 mas in the red and blue 
channels, respectively. The diffraction limit X/D 
is 63 mas and 44 mas for the effective aperture di- 
ameter D ~ 7.5 m set by the coronagraphic mask. 

Figure [T] shows the mosaic of 7 close sub- 
systems found around some primary components 
of wide binaries. Most companions are rather ob- 
vious, well above the detection threshold. The 
PSFs contain also a faint secondary reflex at 13 
pixels (0.24"), with a magnitude difference in the 



■^The oversized spider mask has been taken into account and 
yields a correction factor of 1.079 
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24711 A 



X 



43947A 



43947B 



44804A 



1" 



45734A 



\ 



49520B 



4 

w 

64498A 



ghost 



37332 



Fig. 1. — Mosaic of resolved-component images in the r-channel, labeled by the Hipparcos numbers. Each 
frame is 81 pixels (1.46") across, oriented North down and East to the right. The intensity scaling is negative 
logarithmic. The first 7 frames show all binary components where we detected a close companion. The last 
frame shows the PSF of the unresolved target HIP 37332, with the ghost marked. The faint companions to 
HIP 24711 A and HIP 49520B are marked by arrows. Note that at good AO correction the typical features 
of the NICI PSF become visible. These are particularly the double diffraction spikes and the four spots on 
the first Airy ring, giving it a "boxy" look. These features are caused by the oversized vanes of a spider 
mask placed in the pupil plane to control scatter. 



red channel of about 4.3™. This ghost (illustrated 
m Fig.fflin the case of HIP 37332) is detected in all 
targets, and can be taken advantage of as a fidu- 
cial (see Section 3.4). Its position is the same in 
both channels, whereas the static speckle pattern 
scales in proportion to the wavelength. 

3.3. Measurement of binaries 

To measure the relative position and flux ratio 
of wide binaries more precisely, we fit the scaled 
and shifted image of the primary component (A) 
to the secondary (B). Such fits work very well, 
leaving only small residuals, and produce mutu- 
ally consistent measurements in the two channels. 
The rms scatter between the wide-binary positions 
measured in two NICI channels is 5.5 mas in both 
radial and tangential directions. 

We checked the pixel scale and detector orien- 
tation by comparing the measured positions of 16 
wide pairs (the widest in Table [2] after exclusion of 
HIP 21963, 43947 and 50883, which were not used) 



with the positions listed in the Hipparcos D ouble 



and Muhiple Star Catalogue (|ESA M1997I ). We 



did not find any systematic effect above noise and 
therefore used the nominal pixel scale of 18.0 mas 
for both channels. We found that the detector ori- 
entation of the channels differs by 1.01° ± 0.01°, 
so 1.0° were added to measured angles in the blue 
channel. 

For close pairs where the PSFs overlap, we use a 
variant of iterative blind deconvolution. Good pre- 
liminary estimates of the relative component's po- 
sition and intensity are already available from the 
Gaussian fits of the PSF cores. The image Fourier 
Transform (FT) is divided by the FT of the binary 
with preliminary parameters to get an estimate of 
the PSF. Negative parts of this PSF are set to zero 
and it is replaced by its azimuthal average at ra- 
dial distances over 5 pixels from the center. The 
binary parameters are refined by fitting again the 
image with this synthetic PSF. Then the process is 
repeated, obtaining a second approximation to the 
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PSF and using it again (without azimuthal aver- 
aging) to obtain the final binary parameters. The 
success of this procedure is verified by inspecting 
the resulting PSF, which should be "clean," with 
the least possible trace of the companion. 

3.4. Detection limits 



for our purpose and has been tried be fore with 



good results (e.g. iTokovinin et al. 1120061 ) 




0.2 0.4 0.6 

Radius, arcseconds 

Fig. 2. — Estimated 5a companion detection limit 
(solid line) in the red channel for HIP 37332. The 
simulated companions are over-plotted as asterisks 
(detected) or empty squares (missed). The dotted 
line shows the 2-segment model. 

To determine detection limits, we use mostly 
the data from the red NICI channel where the de- 
tection of low-mass components is deeper because 
of the smaller magnitude difference and larger 
Strehl ratios. The maximum magnitude difference 
of detectable companions Amiim is estimated by 
calculating the rms flux variations (t(p) in circular 
zones of increasing radii p centered on the main 
component. It is assumed that a detectable com- 
panion will have maximum intensity exceeding Scr, 
therefore 



Amiim(/5) < -2.51ogio[5cr(p)//„ 



(1) 



This assumption is, of course, simplistic. We do 
not account for the existence of two simultane- 
ous images in red and blue channels and for the 
good knowledge of the PSF from other companions 
of the same target or from other targets observed 
shortly before or after. Nevertheless, this simple 
strategy of estimating detection limits is sufficient 



The limits computed for HIP 37322 are shown 
in Fig. [21 A bump is produced in the curves at 
p « 13 pixels by the ghost and by fixed speckle 
pattern (see the PSF in Fig. [ij. To confirm the 
validity of ([T|), we performed Monte-Carlo simula- 
tions, adding artificial companions within ±1.5™ 
from the estimated Amn^. Secure detection of the 
ghost (Am — 4.3, p — 0.24") in all images gives 
confidence in our estimation of Amiim. 

The 5cr detection limits Amii„i(p) were fitted 
by two straight lines intersecting at p = 8 pix- 
els (0.144"). The second segment extends from 8 
to 50 pixels (0.9"), but the region between 8 and 
16 pixel radius affected by the ghost and strong 
speckle is avoided in the fit. At p > 50 pixels 
Amiim is assumed to be constant. Only two num- 
bers Am(8) and Am(50) adequately describe the 
actual Amiim (p) curves, with typical rms error of 
only 0.3™ (excluding the zone between 8 and 16 
pixel radii where the 5a method does not properly 
account for the static speckle). 

Detection limits were determined by the above 
procedure for all targets, although the Monte 
Carlo simulation checks were only performed for 
a subset. The two detection parameters Am(8) 
and Am(50) vary little between targets. Their 
median values are 5.02 and 7.79, respectively; the 
quartiles are (4.29, 5.22) and (7.03, 8.06). In the 
following we apply median detection parameters 
to all targets, including the three coronagraphic 
images where the actual limits should be deeper. 

The detection limits found here match other 
similar studies done with "standard" AO. For ex- 
ample, companion detection was complete down 
to AK < 3.9 at 0.14" a nd AK < 7.9 at 0.9" 
in (jTokovinin et al. 2006). Similar r esults were 
obtained bv^Eg genberger et al. ( 2007 ). The deep 



co mpanion survey at Palomar and K eck conducted 
bv iMetchev fc Hillenbrandtl (|2009l) also reached 
AK < 8 at 0.9". Note that typical detection hm- 
its with NICI are much deeper for dedicated planet 
search programs. In these cases long exposure 
times and techniques such as angul ar and spec- 



tral differential imaging ar e applied (jBiller et al. 



20081 : lArtigau et al. Ilioosh . 
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Table 2 

Measured parameters of known and new pairs 



HIP 


Comp 


Red channel (2.272 /nm) 


Blue channel (1.587 ;:im) 


Rem 






e 


P 


Am 


e 


P 


Am 










n 


mag 




n 


mag 




10579 


a,b 


288.58 


6.710 


0.47 


288.50 


6.703 


0.71 




20552 


A,B 


247.35 


5.405 


0.16 


247.55 


5.403 


0.26 


sat? 


21963 


A,B 


92.33 


8.184 


2.35 


92.26 


8.169 


2.64 




23926 


A,B 


259.93 


10.255 


2.35 


259.96 


10.246 


2.66 


sat 


24711 


Aa,Ab 


17.70 


0.670 


5.28 


17.71 


0.663 


5.66 


new 


24711 


B,C 


324.85 


11.135 


6.29 


324.87 


11.132 


6.45 


opt 


27922 


A,B 


19.35 


10.742 


1.55 


19.37 




1.88 




28790 


A,B 


215.28 


5.945 


1.29 


215.18 


5.941 


1.48 


sat 


30158 


A,B 


7.48 


7.096 


1.18 


7.38 


7.093 


1.33 




32644 


A,B 


160.49 


5.029 


0.72 


160.43 


5.039 


0.73 


sat? 


37735 


A,B 


128.43 


6.210 


1.57 


128.49 


6.200 


1.76 




39409 


A,B 


329.64 


5.187 


0.10 


329.67 


5.193 


0.14 




43652 


A,Ca 


128.47 


9.206 


6.16 


128.48 


9.188 


6.41 


opt 


43652 


Ca.Cb 


328.86 


0.124 


0.67 


330.90 


0.126 


0.64 


opt 


43947 


Aa,Ab 


260.79 


0.424 


0.31 


260.78 


0.423 


0.30 


new 


43947 


Ba,Bb 


152.16 


0.289 


1.65 


152.16 


0.289 


1.73 


new 


43947 


Aa.Ba 


202.89 


10.377 


0.82 


202.78 


10.365 


0.90 




44804 


Aa,B 


136.01 


7.175 


1.16 


136.08 


7.167 


1.21 




44804 


Aa,Ab 


313.09 


0.453 


1.84 


313.09 


0.452 


2.02 


new 


45734 


Aa,Ab 


132.69 


0.119 


0.38 


133.33 


0.117 


0.36 




45734 


Aa,B 


193.63 


9.077 


0.71 


193.65 


9.076 


1.23 




45940 


A,B 


69.28 


6.713 


1.72 


69.29 


6.705 


1.98 




47836 


A,C 


328.29 


9.809 


5.31 


328.36 


9.797 


5.77 


opt 


49520 


A,B 


326.95 


9.552 


0.29 


326.95 


9.543 


0.44 




49520 


Ba,Bb 


44.21 


0.213 


3.61 


46.69 


0.210 


3.40 


new 


50883 


A.B 


192.65 


6.839 


2.33 


192.57 


6.830 


2.57 




50883 


A,C 


299.87 


7.379 


6.23 


299.82 


7.363 


6.55 


opt 


55288 


A,B 


253.94 


9.462 


0.54 


253.96 


9.456 


0.84 




58241 


A,C 


118.32 


4.015 


6.38 


118.09 


4.008 


6.52 




64498 


Aa,B 


302.20 


9.464 


1.59 


302.24 


9.458 


1.91 




64498 


Aa,Ab 


296.55 


0.372 


1.76 


296.31 


0.370 


2.01 


new 


65176 


A,B 


96.89 


5.168 


0.61 


96.95 


5.164 


0.62 
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4. Results 

4.1. Companion data 

We have found 6 new close companions, one 
previously known companion, and several faint op- 
tical companions (Fig.[T]). Table[2]lists the relative 
positions and magnitude differences in known and 
new pairs. Each pair is identified by the HIP num- 
ber and component designations. The measured 
position angles in degrees, separations in arcsec- 
onds, and magnitude differences are listed sepa- 
rately for the two NICI channels. 

The large rms difference between Hipparcos and 
our measurements of the wide pairs (130 mas in 
p and 1.05° in 9) is mostly caused by their mo- 
tion between the two epochs. The comparison 
with the relative positions of the same binaries 
deduced from the 2MASS Point Source Catalog 



(jCutri et al. II2003I ) has confirmed that no system- 



atic corrections are needed and has shown a scatter 
of similar magnitude. 

Remarks in the last column of Table [2] indicate 
cases when the image was saturated (the relative 
photometry may be compromised) , newly detected 
close companions, and likely optical companions. 
These faint companions with separations of a few 
arcseconds would be dynamically unstable within 
wide binaries, should they be physical!. 

4.2. Notes on individual systems 

HIP 28790 has a common-proper- motion com- 
panion C at 196". 

HIP 36165 A is suspected to be a spectro- 
scopic binary (SB). 

HIP 43947 is unexpectedly discovered to be 
a resolved quadruple system. Both components 
are located ~ 2™ below the Main Sequence in the 
{Mv, V — K) CMD, possibly because the Hippar- 
cos parallax measurement was biased by the mo- 
tions in the sub-systems. The position and Am of 
the wide pair Aa,Ba is derived from the Gaussian 
fits, therefore it may be less accurate. 

HIP 45734: This is a young pre-Main- 
sequence system and X-ray source RX J0919.4- 



Dynamical stability of a triple system is possible when the 
ratio of the inner semi-major axis to the outer p eriastron 
distance is larger than ^3 fe.e. lHarrington"1ll972l ). As only 
projected separations are known for our triples, assessment 
of their dynamical stability is approximate. 



7738. It projects on the Chamaeleon star-forming 
region, but is located much closer than the corre- 
sponding association. ICovino et al. ] (Il99i1) found 
from a single spectrum that component B (south- 
ern) is a double-lined SB, with a large RV d iffer- 



ence from A (see also Desidera et al. 2006[) . A 



itself is a close visual b i nary KOH83 resolved in 
March 1996 bv iKohler I (|200l[ ). He measured the 
position (173.9°, 0.109") and AK = 0.45. This 
is the only previously known resolved sub-system 
among the observed stars. Follow-up observations 
will soon lead to the visual orbit and mass mea- 
surement, adding another empirical point to the 
evolutionary tracks of young stars. 

HIP 49520: The faint companion Bb is close 
to the detection limit (see Fig. [Ij. Yet, it is se- 
curely detected in both channels and stands up 
clearly in residuals when fitting the image of B 
with the image of A. 

HIP 50638/36: The wide pair has been ob- 
served since 1835; it is definitely physical. Compo- 
nent B is variable. For this reason, probably, it is 
located on the CMD some 6™ below the MS, while 
A is on the MS with the parallax t thip = 83 mas 
listed in Table [TJ According to Ivan Leeuwen 



(|2007l ). the parallaxes and proper motions of the 
two components are discordant, but have large er- 
rors. Future RV measurements can settle the con- 
troversy. 

HIP 59021: Component A is a suspected SB 
and also an astrometric binary (acceleration so- 
lution in Hipparcos). Non-detection of the sub- 
system with NICI means that its period is likely 
less than lOy, or it is a white dwarf. 

4.3. Estimation of the mass ratios 

iC-band absolute magnitudes Mk of all com- 
panions were computed based on the photom- 
etry from 2MASS (see Table [IJ and trigono- 
metric parallaxes tthip from Hipparcos. For re- 
solved pairs, we assume that the Am in the 
red NICI channel correspond to magnitude dif- 
ferences in the if-band. The masses were es- 
ti mated from Mk using the empirical relations 
of iHenrv fc McCarthvl (|l993l ). We also tried a 
quadratic app roximation of relation between mass 
and Mk from lCirardi et al. I (l2000l) The relative 



difference with iHenrv fc McCarthv I ( 1993 ) is less 
than 10% for masses above 0.1 M©. We checked 
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Separation, AU 

Fig. 3. — Mass ratios q — M2/M1 of 7 the re- 
solved sub-systems versus their projected separa- 
tions. The shaded region indicates incomplete de- 
tections (light grey corresponds to detection prob- 
abilities from 0.5 to 0.9, dark grey from 0.1 to 0.5). 



that the influence of using cither of these rela- 
tions on the final result is not critical. Table [3] 
lists the estimated masses of the companions in 
the resolved sub-systems, their projected separa- 
tions r = /o/tthip and order-of-magnitude period 
estimates P* = r^/'^{Mi + Mi)-!/^. 

Figure [3] shows the mass ratios q = M2/M1 
in the resolved inner sub-systems as a function 
of projected separation. The median parallax 
of observed targets is 19 mas, so the upper limit 
of the separation range, 100 AU, corresponds to 
p w 1.9". Maximum separations of dynamically 
stable sub-systems should not exceed 100 AU in 
most cases because the projected separations of 
the wide binaries range from 150 to 1000 AU (Ta- 



Table 3: Estimated parameters of the 7 observed 



HIP 


Comp 


Ml 


M2 


r 








Me 




AU 


yr 


24711 


Aa,Ab 


1.04 


0.15 


43.6 


260 


43947 


Aa,Ab 


0.64 


0.59 


15.9 


60 


43947 


Ba,Bb 


0.55 


0.24 


10.9 


40 


44804 


Aa,Ab 


0.95 


0.61 


29.3 


130 


45734 


Aa,Ab 


0.88 


0.81 


6.8 


14 


49520 


Ba,Bb 


0.92 


0.26 


12.6 


40 


64498 


Aa,Ab 


1.09 


0.71 


20.6 


70 



ble [T]). Indeed, all detected sub-systems have 
r < 44 AU. The lower limit of 5AU corresponds 
to p ~ 0.1", well above the NICI resolution limit. 

The detection limits in the red channel were 
converted into mass ratios and averaged over the 
61 observed components, producing a smooth de- 
tection probability p^^t [f, q) ■ Contours of this 
function at 0.1, 0.5, and 0.9 levels are over-plotted 
in Fig. [21 We see that the chosen region of the 
(r, q) parameter space is well covered with NICI. 

4.4. Statistical analysis 

The purpose of this study is to reach statisti- 
cal conclusions about the frequency and the prop- 
erties of sub-systems in solar-type wide binaries. 
Considering the small number of components, we 
keep this analysis at a basic level, and assume that 
the mass ratios in secondary sub-systems have a 
power-law distribution independent of separation: 



/(g) 1)/. 



(2) 



This model has only two free parameters, the total 
companion frequency e and the power-law index /3. 
Neglecting incomplete detection, the companion 
count gives e 7/61 = 0.12 ± 0.04. 

The distribution of the companion projected 
separations in the considered range from 5 to 
100 AU matters only because the detection limits 
depend on r. We assume that /(logr) = const. 

The input data are the (r, q) values for K = 7 
sub-systems, the total number of surveyed targets 
iV = 61, and the array of detection probabilities 
Pdot('', ?)• We find the estimates of parameters 
(e, f3) by the Maximum Likelih ood (ML) method 
(see e.g. iTokovinin et al. Iliooei Apendix B). The 
likelihood function £ is a product of the proba- 
bilities to obtain observations for each component 
(non-detections ior N — K systems or detections 
of K subsystems with mass ratios qk)- 

The minimum of 5 = — 21n£ is sought. This 
function depends on the input data and the model 
parameters as 

K 

S = 2Nfo - 2^1n[/(gfc)pdet(rfc,(Zfc)]. (3) 

fc=i 

The summation over k is done for the 7 detected 
companions. The companion probability /o equals 
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the product /(r, q)pdcti'r, q) averaged over the con- 
sidered part of the parameter space {r,q). If the 
detection probability pdot equals one everywhere, 
the result is fa ~ e, otherwise /o < £• 



-5 0.5 - 




-0.5 



0.1 0.2 
Companion frootion 



0.3 



Fig. 4. — Contours of the ^-function in the e, f3 
space corresponding to the 68% and 90% confi- 
dence areas. The minimum is marked with a cross. 



The contours of S in the parameter space define 
confidence limits: AS = 1 corresponds to the 68% 
interval ("la"), AS* = 2.71 to 90% and AS = 4 
to 95%, in direct analogy with the Gaussian prob- 
ability distribution. We verified the ML method 
on a simulated data set. The input parameters 
were recovered and the confidence intervals were 
as expected. 

Minimization of ([3]) leads to (e, P) = (0.12, 0.18). 
If we ignore the detection limits and set pdot = 1, 
the result changes little: (e, /3) = (0.11, 0.40). The 
companion fraction 0.12 is not different from its 
naive estimate. As to the power index /3, the 
ML points to a uniform distribution in q as being 
most likely, although the uncertainty is large. The 
confidence areas of the parameters are shown in 
Fig. m Considering that the contours are not very 
elongated, we can determine the 68% and 90% 
confidence intervals for each parameter by a simple 



Table 4: Confidence intervals 



Interval 


e 


/3 


Min. S 


0.12 


0.18 


68% 


0.08.. .0.17 


-0.30. ..0.80 


90% 


0.06.. .0.21 


-0.54.. .1.28 



cross-section through the minimum (TableU]). The 
data are marginally compatible with /3 ~ —0.5, in 
which case the fraction of sub-systems may ex- 
ceed 0.2 (see the lower-right extension of the 90% 
contour in Fig. 



5.1 



Discussion 
Completeness of the existing catalog 



The fact that, in a sample of only 33 sys- 
tems, six of the seven observed companions are 
new, shows that current knowledge of multiplicity 
among nearby solar-type systems is very incom- 
plete, and that near-IR AO imaging is a powerful 
technique for detecting new companions. 



40 F 



30 



20 



10 



* 



^ 



L_ 
0.0 



* 



0.4 0.6 
Moss ratio 



0.8 



1 .0 



Fig. 5. — Histograms of the mass-ratios in inner 
systems of triple stars. Solid lines correspond to 
the 123 close pairs and dashed lines to the 61 wide 
pairs. The vertical bars show ±^/N statistical 
errors for each bin of width 0.2. All points are 
slightly displaced horizontally for clarity. 

We plot in Fig. [5] the mass-ratio distribution 
in the inner sub-systems of known triple stars. 
The data are extracted from the upda ted version 
of th e Muhiple-Star Catalog (MSC, iTokovinin I 
19971 ). where component masses are estimated by 
a variety of methods. We selected 220 systems 
within 200 pc from the Sun with masses of pri- 
mary components from 0.5 to 1.5 Mq. They were 
sub-divided by the separation in the inner sub- 
systems in two groups, 123 close (semi-major axis 
< lAU) and 61 wide (between 5 and 100 AU). 
The remaining 36 pairs have even wider separa- 
tions. The mass-ratio distributions in these two 
groups arc markedly different, but our un-biased 
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survey results in a flat mass-ratio distribution, 
within errors. Therefore, the difference seen in 
Fig. [5] is caused by observational selection. Given 
an approximately equal number of cfose and wide 
sub-systems with large mass ratios (two last bins 
in Fig. [5]) and assuming that the mass-ratio dis- 
tributions in these two groups are indeed simi- 
lar, we expect that the number of close and wide 
sub-systems (with separations below lAU and 
between 5 and 100 AU, respectively) are similar, 
too. The completeness of the MSG for wide sub- 
systems is thus about 1/2. 

5.2. Companion fraction 

Is the presence of a wide binary companion in- 
fluencing the frequency of inner sub-systems? To 
answer this question, we compare our results on 
triples with surveys of solar-type binaries. Unfor- 
tunately, binary-star statistics are usually derived 
regardless of higher-order multiplicity, so a cleaner 
comparison must await new study of a large unbi- 
ased sample. 

The separation range from 5 to 100 AU explored 
here corresponds to orbital periods between 10 y 
and 1 000 y. According to iDuauennov fc Mavor 



(|1991I ). about 20% of solar- type dwarfs have com- 
panions with such periods. In the actual, nar- 
rower range of projected separations from 6.8 AU 
to 44 AU, 12% of solar-type dwarfs have com- 
panions. Therefore, the frequency of sub-systems 
around components of wide binaries is similar to 
the frequency of binary companions to single stars 
in the same separation range. 

The frequency of sub-systems with periods be- 
low 3jHii_j2^ual_jnult^^ has b een determined 
by iTokovinin fc Smekhov (2002) to be between 
11% and 18%. Again, it turned out to equal 
within observational errors the companion fre- 
quency among single dwarfs of similar masses in 
the solar neighb o rhood and in some open clusters. 
Desidera et al. ( 20061) surveyed RVs of 56 wide 
visual pairs and found that 0.135 ± 0.05 fraction 
of companions have spectroscopic sub-systems. 

Therefore, one might surmise that the presence 
of a wide companion does not affect the formation 
of inner, closer sub-systems in a wide range of sep- 
arations. This question is however far from having 
been settled for both stellar and planetary sub- 
systems (see the discussion in lEggenberger et al. 



20071 ) . We know that almost all close binaries with 
periods below 3 d do have additional outer com- 
panions, although the mass ratios of close binaries 
with and without tertiary companions ha ve iden- 



tical distributions ( Tokovinin et al. 1120061 ). 



5.3. 



Is the mass-ratio distribution univer- 
sal? 



Metchev fc Hillenbrandt I (|2009l) advocate the 



idea of a companion mass function with a uni- 
versal power law over a wide range of separations 
and masses. They fit the mass-ratio distribution 
in solar-mass visual binaries to a power law with 
/? — — 0.39±0.36. Visual companions to stars more 
massive than the Sun follow a power law f{q) oc 
with /? from —0.3 t o —0.5, as established in AO 



ima ging surveys bv IShatskv fc Tokovinin I (|2002l ) 
and Kouwenhoven et al. ()2005 ). Gertainly, these 
distributions do not correspond to the most sim- 
plistic case of rand om companion pairin g as d e- 
scribed in detail bv iKouwenhoven et al. ] (|2009[ ). 

Our survey of triples points to a flat mass- 
ratio distribution with /? ^ 0. The errors are 
large, so the difference with visual binaries is not 
yet significant. The mass-ratio distribution in 
closer, spectroscopic binaries is approximate ly flat 



(|Mazeh et al. II2003I : Itfalbwachs et al. II2003I ). 



ilar to the mass-ratio distribution in close sub- 
systems of triple stars (Fig. [5|) . If we confirm on a 
larger sample that sub-systems in triples do differ 
from pure binaries in their mass-ratio distribution, 
its universality will be seriously challenged. 

5.4. Brown-dwarf desert 

The mass function of companions to binary and 
multiple stars is radically different from the Initial 
Mass Function. Therefore, the fraction of low- 
mass companions is much smaller than the frac- 
tion of low-mass stars in the field. This fact, 
known as brown dwarf (BP) desert, i s well es- 



tablished ( Grether fc Lineweaver 2006[ ). Adopt 



ing a flat mass-ratio distribution for our sample, 
also typical for spectroscopic binaries, we estimate 
that only 7% of companions to a solar-mass star 
will be in the BD regime. Gonsidering the com- 
panion frequency of 0.12 in the studied separation 
range, one BD companion per ^100 primaries is 
expected. Our result thus suggests that the BD 
desert in multiple systems extends to separations 
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of 100 AU. 

The BD desert and the tendency to equal com- 
panion masses are explained in the current star- 
formation scenario by continuing accretion onto a 
binary. Whenever a low-mass companion forms 
by fragmentation of a massive disk, continuing ac- 
creti on usually brings its mass into stellar regime 
(e.g. iKratter et al. 2010l ). Companions of sub- 
stellar masses should be intrinsically rare. At the 
same time the orbital drag (inward migration) re- 
duces the separation, producing close sub-systems. 
It appears that about half of the companions set- 
tle in close orbits with separation < 1 AU, the rest 
of them remain in wider orbits. 

6. Conclusions 

We surveyed 33 nearby wide binaries with solar- 
type primaries and found 7 resolved sub-systems, 
most of them previously unknown. We derive the 
fraction of sub-systems with projected separations 
from 5 to 100 AU to be 0.12 ± 0.04. The mass- 
ratio distribution in these sub-systems appears to 
be flat, to within a large statistical error. 

The sample can be increased 6 times if we con- 
tinue observations with NICI and extend them to 
the Northern sky using some other AO facility. 
Such investment of 5 nights on 8-m telescopes will 
enable to reduce the size of the confidence inter- 
vals in Fig. lUby 2.5 times. Then, it will become 
possible to distinguish between the negative power 
index /? ~ —0.3 typical for resolved binaries and 
flat mass-ratio distribution /? ~ characteristic of 
close binaries. 

The joint fraction of both visual and spectro- 
scopic sub-systems in wide solar-type binaries is 
about 0.25. In consequence, the majority of the 
components are single, and provide an accomo- 
dating environment to host planetary systems. 

We thank M. Chun and F. Rigaut for help- 
ing us diagnosing a NlCl AO misbehaviour, and 
F. Rantakyro for his assistance in the preparation 
and execution of the observations. We are thank- 
ful for discussions on the manuscript to W. Brand- 
ner and N. Huelamo. This work used the NASA's 
Astrophysics Data System, data products from the 
2MASS funded by the NASA and the NSF, and 
the SIMBAD database maintained by the Univer- 
sity of Strasbourg, France. The comments of the 



anonymous referee helped us to improve the pre- 
sentation. 
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